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Complexiton solutions (or complexitons for short) are exact solutions newly introduced 
to integrable equations. Starting with the solution classification for a linear differential 
equation, the Korteweg-de Vries equation and the Toda lattice equation are considered 
as examples to exhibit complexiton structures of nonlinear integrable equations. The 
crucial step in the solution process is to apply the Wronskian and Casoratian techniques 
for Hirota's bilinear equations. Correspondence between complexitons of the Korteweg-de 
Vries equation and complexitons of the Toda lattice equation is provided. 

1. INTRODUCTION 

Differential equations or differential- difference equations can describe various motions 
in nature. It is important to study their integrable properties, and more important to 
tangibly determine their exact solutions. Soliton theory is one of significant developments 
along this direction. The theory tells us that there exist soliton solutions to many in- 
tegrable equations, both continuous and discrete. More generally, negatons (generalized 
solitons) can be explicitly presented (see [CI El, for example, for the Korteweg-de Vries 
(KdV) equation and the Toda lattice equation). Similarly, there exist positons 
which is another important achievement in soliton theory. What could we say about exact 
solutions further? This report is aiming at discussing this question. Specifically, we will 
explore what other solutions to integrable equations can exist, and show that so-called 
complexiton solutions [|7j are one of new exact solutions. 

Let us first observe an example of linear ordinary differential equations: 

ay" + by' + cy = 0, y' — dy/dx, y" = d 2 y/dx 2 , a,b,c = real consts., (1) 

to recall the solution classification of linear differential equations. Its characteristic equa- 
tion is a quadratic equation 



am 2 + bm + c = 



and the quadratic formula gives its two roots: 



-b ± \Jb 2 - Aac 
mi < 2 = 2a • 
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These two roots completely determine the general solution y = y(x) of (0). We now 
describe all solution situations of ((H), with pointing out the corresponding solutions in 
soliton theory. 

• Real roots: 

— Distinct roots ni\ ^ m 2 : The general solution of (0) is then 

y (x) = c ie mix + c 2 e m2X , ci, c 2 G R. 

This corresponds to so-called solitons and negatons in soliton theory. 

— Repeated roots mi = m 2 : The general solution of ([TJ is now 

y [x) = Cl e mix + c 2 xe mix , ci, c 2 G R. 

This corresponds to negatons of higher order when mi = m 2 and rational 
solutions when mi = m 2 = in soliton theory. 

• Complex roots: 

— Purely imaginary roots mi 2 = ±/3y/—l: The general solution of (0) is then 

y(x) = c\ sin(/3x) + c 2 cos(/5x), Ci, c 2 G R. 

This corresponds to so-called positons in soliton theory. 

— Not purely imaginary roots m ii2 = a ± (3yf—l\ The general solution of (JTJ) is 
now 

y(x) = cie ax sm(/3x) + c 2 e ax cos(ftx), ci, c 2 G R. 

This corresponds to so-called complexitons which we are going to discuss. This 
solution can also boil down to periodic solutions of positon type if a — > 0, and 
exponential function solutions of negaton type if (3 — > 0. 

It is always possible to classify exact solutions of constant-coefficient, linear ordinary 
differential equations of any order. A great review of solution classifications was given in 
Ince's book [|H], where the solutions of the second-order linear ordinary differential equa- 
tions were classified in terms of hypergeometric functions, Riemann P-functions, etc. The 
problem of solution classifications becomes very difficult for nonlinear differential equa- 
tions. Galois differential theory may help in handling differential equations of polynomial 
type. We will only concentrate on integrable equations, for which we can start from their 
nice mathematical properties such as symmetries and adjoint symmetries [151 110j. 

What is a complexiton solution? The general notion of complexitons could be 
characterized by the following two criteria: 

• Complexitons involve two kinds of transcendental functions: trigonometric and ex- 
ponential functions. 

• Complexitons correspond to complex eigenvalues of associated characteristic prob- 
lems. 
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This report aims at constructively contributing to the theory of complexiton solutions to 
nonlinear integrable equations, and the KdV equation and the Toda lattice equation will 
be taken as two illustrative examples. The proposed idea of constructing complexitons 
through special determinants, for example, the Wronskian and Casorati determinants, 
will also work for other integrable equations. In particular, super-complexitons can be 
similarly generated for super integrable equations, which will extend the theory of super- 
solitons [QUQ2]. 



2. KORTEWEG-DE VRIES EQUATION 

Let us first consider the KdV equation 

u t + u xxx - Quu x = 0, 



where u+ = %, u r = and u rrr = It is known that under the transformation 



u = -2(ln/) a 



2(/ fxx fx) 

T 2 : 



(2) 



(3) 



the KdV equation (j2J) is transformed into the bilinear equation 
(D x D t + D 4 x )f-f = 0, 
that is, 

fxtf ft fx fxxxxf ^fxxxfx ~t" ^ fxx ^' 

where D x and D t are Hirota's operators: 

oo 1 

f( x + M + k ) g ( x _ K t _ k ) = —(//,/>// • g)tik\ 

I. J. 



(4) 



i,j=0 



or more directly, 



x'=x,t'=t 



A powerful method of solutions for integrable bilinear equations is the Wronskian tech- 
nique [Ej- To construct solutions, we use the Wronskian determinant 



f = w(<i> 1 ,(h 



)■■- 



,(0) 
N 



,(AT-1) 

; v 



(5) 



where 



d 1 * (pi / dxi , j > 0. The resulting solutions are called Wronskian solutions. 



The Wronskian technique requires 



1 < i < JV, 
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and thus all involved functions 0j, 1 < % < N, are eigenfunctions of the Lax pair of 
the KdV equation associated with zero potential. Actually, the Wronskian solution can 
be generated from the Darboux transformation of the KdV equation starting with zero 
solution. The above system generates the eigenfunctions needed in forming Wronskian 
solutions: 



2i; 



a" 



4>i = Cux + G 

4>i = C u cosh(r)iX - irjft) +C 2i smh(r)iX - irjft), r]i 
<pi = C u sin(?7jX + 47/ft) + C 2 i cosfax + rjj = 

when Aj is zero, negative and positive, respectively. Here Cu and C 2 i are arbitrary real 
constants. 

More general Wronskian solutions can be constructed under a broader set of sufficient 
conditions [I14j: 



= £ A 

3=1 



H,t 



-Hi 



1 < i < N, 



where the coefficient matrix A := (Ay) is an arbitrary constant lower-triangular matrix. 
Very recently, an essential generalization to the above sufficient conditions is presented in 



N 
3=1 



Mi 



1< i < N, 



(6) 



where the coefficient matrix A = (Ay) is an arbitrary real constant matrix, not being lower- 
triangular any more. Once a set of eigenfunctions is obtained, the Wronskian solutions 
to the KdV equation (j2J) is given by 



u 



-2fl£lnW(0 



i,<?2, ••-,9n) 



(7) 



It is easy to see that linear transformations of eigenfunctions do not generate new 
Wronskian solutions. Therefore, we only need to consider the Jordan form of the coefficient 
matrix A. A real matrix can have and only have two types of Jordan blocks in the real 
field: 



Type 1: 



A, 
1 







A; 



Type 2: 







1 A,; 



Ai 







' 1 


" 




, h = 






_ Pi a { 
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where Aj, ai and > are real constants. The first type of blocks only has the real 
eigenvalue Aj with algebraic multiplicity hi, but the second type of blocks has the complex 
eigenvalues Xf = ai ± /3j\/^T with algebraic multiplicity Zj. 

Note that an eigenvalue of the coefficient matrix A = (Ay) is also an eigenvalue of 
the Schrodinger operator — ^ + u with zero potential u = 0. Therefore, according to 
the types of eigenvalues of the coefficient matrix A, we can have four Wronskian solution 
situations [I15j: 



Rational solutions: 


Type 


1 


with 


A, = 0, 


Negaton solutions: 


Type 


1 


with 


A, < 0, 


Positon solutions: 


Type 


1 


with 


A, > 0, 


Complexiton solutions: 


Type 


2 


with 


complex eigenvalues. 



Complexiton Solutions of Zero Order: Assume that 



4>il,xx 






j Ai — 


O-i ~Pi 


1 


<M,t 


= -4 


YiX,xxx 


<fii2,xx 




. ^ 2 - 




Pi &i 




. ^ 2 >* . 




4 > i2,xxx 



(8) 



where a^ and /3j > are arbitrary real constants. Such eigenf unctions can be easily 
explicitly presented [ [7j, but it is not easy to get the eigenfunctions associated with 
higher-order Jordan blocks of type 2 (see [E3 for more information). An n-complexiton 
solution of order ■ ■ ■ , l n ) is defined by 
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-2d 2 x ln W(<f> u , 012, ^11,^12; 



• " • ) <Pnl, 



^■■■,d l «J nl ,d l 



>n2) 



which corresponds to Jordan blocks of type 2. An n-complexiton of order (0, ■ 
complexiton for short) is 

U = -2d 2 x ln W((pii, 012, ■ • ■ , 0nl, <pn2)- 

In particular, a 1-complexiton [|7j is given by 

u = -2^]nW(<hi,<hi) 

-4pf[l + cos(25 1 (x - Pit) + 2k x ) cosh(2Ai(a; + a t t) + 2 7l )] 

[Ai sin(25i(x - Pit) + 2«i) + <5i sinh(2Ai(x + out) + 271)] 2 
AaiPi sin(25i(x - Pit) + 2«i) sinh(2Ai(x + out) + 271) 
[Ai sin(25i(x - Pit) + 2«i) + Si sinh(2Ai(:r + ait) + 271)] 2 ' 

where oti, Pi > 0, Ki and 71 are arbitrary constants, and Ai, Si, di, and Pi are 



(9) 
,0) (n- 



Ai 



'a{ + P{ — ai 



Si 




n, =4 A /a 2 + /5 2 + 8ai, ft = 4y«f + Pi - %a x . 
The special case of at\ = leads to the following solution 

8/3i + 8/3i cos( v / 2A x - 4/3 lv ^t + 2«i) cosh(v^x + Ap l ^Wit + 2 7 i) 



u 



[sin(v^x - 4PiV2~Pi~ t + 2«i) + smh{^/2p^x + Ap l ^2p[t + 271)]' 
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This solution is associated with purely imaginary eigenvalues of the Schrodinger spectral 
problem with zero potential. Moreover, the 1-complexiton above contains the breather- 
like or spike-like solution presented in [llfij: 

8{(a 2 — b 2 )(b/a) cosz/sinh(?7 + p) + 26 2 [1 + sinz/cosh(?7 
[cos v — (b/a) sinh(?7 + p)] 2 

where v = —2bx + 8(3a 2 6 — b 3 )t, T) = —2ax + 8(a 3 — 3ab 2 )t and p = ln(6/a). This is 
obtained if we choose that 

7T 1.6. , 

4 2 a 

where a and 6 are arbitrary real constants. Two special cases of the 1-complexiton solution 
are depicted in Figure 1. 




Figure 1. 1-complexiton - a,\ = 0, /3i = 1 (left) 1-complexiton - a\ = 1, j3\ = 1 (right) 



Complexiton solutions of higher order are complicated but special ones can be presented 
by taking derivatives of eigenfunctions with respect to the involved parameters [[7J IT5] . 
Complexiton solutions are singular and they are not travelling waves. It is clear that 
the 1-complexiton above is not a travelling wave since (3\ 7^ 0. Moreover, based on the 
Painleve property of the KdV equation, the singularities of complexitons are all poles of 
second order with respect to x. This is obvious for the 1-complexiton above, since u has no 
pole singularity (the numerator of u will be zero) if the function Ai sin £1 + 5i sinh £ 2 is zero 
and its spatial derivative 2A 1 <5 2 (cos£i + cosh£ 2 ) is also zero, where £1 = 25\{x — j3\t) + 2ki 
and £ 2 = 2A 1 (a: + a x t) + 271. 

3. TODA LATTICE EQUATION 

Let us second consider the Toda lattice equation 
a n = a n {b n -i - b n ), b n = a n - a n+1 , (10) 
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where a r 



^ and b n = Under the transformation 



a " = 1 + 3^ lo g r ™ = " 



& « = Til log — 
r n+l 



d * 2 ~ r n 

the Toda lattice equation (fTUj) becomes 

[D t 2 -4sinh 2 (^)]r n .r ft = 0, 

where D t and _D n are Hirota's operators. That is, 



T n T n+l 



'ID 



Tn+l T n-l + r n 



0. 



(12) 



The Casoratian technique [ Ej is one of the ways to solve the bilinear Toda lattice 
equation ([12]). The corresponding solutions to the bilinear Toda lattice equation (jT2j) are 
determined by the Casorati determinant: 



&i(n + JV- 1) 
&a(n + JV- 1) 



0i (n) 0i(n + l) 
2 (n) 02(^ + 1) 

r n = Cas(0i,0 2 ,---,0A r ) : 

07V (n) <p N (n + l 
and the Casoratian technique requires [ ITS] 

JV 

0i(n + 1) + 0i(n - 1) = ^ Ajj^j(n), <9 t 0j(n) = 0j(n + 5), 1 < i < iV, 



(13) 



(14) 



where 5 = ±1 and \j are arbitrary real constants. Similarly, we only need to consider the 
Jordan form for the coefficient matrix A = (A^-). Assume that two types of Jordan blocks 
of the coefficient matrix A = (Xij) are specified as in the case of the KdV equation in 
the last section. Then, based on the solution structures of eigenfunctions associated with 
eigenvalues of different types, we can similarly have four Casoratian solution situations: 



Rational solutions: 


Type 


1 


with 


Ai = ±2, 




Negaton solutions: 


Type 


1 


with 


Xi < -2 or 


Xi > 2, 


Positon solutions: 


Type 


1 


with 


-2 < Aj < 


2, 


Complexiton solutions: 


Type 


2 


with 


complex eij 


^envalues. 



Case ki = 1 of Type 1: We have 

0j(n + 1) + 4>i(n - 1) = Aj0j(n), d t (f>i(n) = 4>i(n + 5), 

where 5 = ±1 and Aj = consts. Its general eigenfunctions are: 

fa = C u (Sn + et)e n e £t + C 2i e n e Et , A, = 2e, 

= Ciid tcoskl cos(kin + 8t sin ki) + C 2l e tcoskl sin(kin + 5t sin ki), Aj = 2 cos ki 



(15) 



Oijt e -+- 02j£ e 



A, = 2e cosh ki, 
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where e = ±1, | cosfcj| < 1 for the second type of eigenf unctions, cosh/cj > 1 for the third 
type of eigenf unctions, and Cu and Ci% are arbitrary real constants. These three types of 
eigenfunctions lead to rational solutions (see, for example, [HE]), negatons and positons 
(see, for example, [121), respectively. 

Case li — 1 of Type 2: To construct complexitons, we solve 



(16) 



<f>i(n + 1) + <Px(n - 1) = «0i(n) - /3<j) 2 (n), ( d t 4> x {n) = (f>i(n + 5), 
<f> 2 (n + 1) + 4>2(n - 1) = /90i(n) + a0 2 (n), { <9 t 2 (n) = 2 (n + 5), 

where 5 = ±1 and a, (3 — consts. Assume that 

a + PV~L = A = 2 cosh fc = e k + e" fc , k = a + b\T-\ , 

where a and b are real, and then 

a = 2 cosh a cos b, (3 = 2 sinh a sin 6. 

View (JTBJ) as a compact equation like (fT5|l where A = a + an d then using the third 

type of eigenfunctions above, we can generate the corresponding eigenfunctions of ()16|) : 

^(n) = C e« a + 4e5 ' lcos ' 5 '' cos(n6 + £e 5 ° sin 56) + c fe- na+ * e ~ 5acos<5b cos(nfe + te~ 6a sin 56), 

^ a (n) = ce na+e5acosSb sm(nb + te Sa sm5b) - rf e - na+ * e_aacos5fc sm(nb + te~ Sa sm5b). ^ 

This set of eigenfunctions yields the r-function of 1-complexiton: 

T n = Cas(0i(n),0 2 H) 

= 2cd e 2i cosh Sa cos 5b sin (2n6 + b + 2t cosh 5a sin 5b) sinh a 

where 5 = ±1 and a, b,c,d— consts. In particular, c = ±d leads to the r-function: 

r n = 2c 2 e 2tcosh5acos5b smh(2na + a + 2tsmh5acos5b)smb 

±2c 2 e 2t coshSacosSb S in{2nb + b + 2t cosh 5a sin 5b) sinh a . 

The resulting 1-complexiton solution of the Toda lattice equation (jTOj) is given by (fTTJ) 
where r n is presented in f)18jl . Two special cases with 5 = c = d= lof this 1-complexiton 
solution are depicted in Figures 2 and 3. Complexiton solutions of higher order are com- 
plicated but special ones can be constructed by computing derivatives of eigenfunctions 
with respect to the involved parameters [I19j. 

4. CONCLUDING REMARKS 

There is a correspondence between the characteristic linear problems of the KdV equa- 
tion and the Toda lattice equation: 

d 2 x & 0(n + l) + 0(n-l)-20(n), 
-dl<f) = X(f) 4» 4>(n + 1) + <f>(n - 1) = \<f>(n), 
-d 2 x (j) = \(j) -<j>(n + 1) -</>(n- 1) + 2<j)(n) = (-A + 2)0(n), 
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Figure 2. 1-complexiton a n - 2a = b = 2 (left) 1-complexiton b n - 2a = b = 2 (right) 




Figure 3. 1-complexiton a n - a = — b = 1 (left) 1-complexiton b n - a = — b = 1 (right) 

and a correspondence between the eigenvalues of the characteristic linear problems: 

A = 0, > 0, < and A = complex |A| = 2, < 2, > 2 and A = complex, respectively. 

We can also observe the correspondence between the KdV equation and its characteristic 
linear problem — <9^0 = A0. The characteristic equation of —d^cf) = \<j) is 

-m 2 - A = 0. 

Therefore, for example, we see that 

A is positive <^ m is purely imaginary, 

A is complex <^ m is complex but not purely imaginary. 
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Together with other obvious correspondences, this implies that 
A = KdV rational solutions, 
A < <^> KdV solitons and negatons, 
A > <^> KdV positons, 
A = complex <^ KdV complexitons. 

It then follows that the obtained complexiton solutions to the KdV equation and the 
Toda lattice equation satisfy the two criteria for complexiton solutions stated in the 
introduction, indeed. Because of the characteristic in the second criterion, complexitons 
may play a role similar to the one that the imaginary unit plays in physics. 

We also mention that the Wronskian determinants and the Casorati determinants give 
many other solutions if different types of eigenvalues are allowed. Such solutions are 
called interaction solutions among determinant solutions of different kinds. For higher 
dimensional integrable equations, the solution situations are much more diverse [ I2H1 E] 
and the problem of classifying solutions is extremely difficult [122]- It is hoped that the 
study of complexitons could further assist in understanding, identifying and classifying 
nonlinear integrable differential equations and their exact solutions. 
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